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Abstract 

Density perturbations in the flat {K = 0) Robertson-Walker universe with radiation (p = e/3) and 
positive cosmological constant (A > 0) are investigated. The phenomenon of anomalous dispersion of 
acoustic waves on A is discussed. 
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1 Introduction 


The image the perturbations create on the last scattering surface depends on how they propagate in the radia- 
tional era. In the flat radiation dominated universe with vanishing cosmological constant the density perturba¬ 
tions form sound waves and propagate with the constant sound velocity v = l/-\/3. They obey linear dispersion 
relation lj = k/y/3, therefore, their phase and group velocities are identical. Their propagation does not depend 
on the wave number, therefore, no critical scale appear and no gravitationally bound structures form. The 
above is the subject of theorem proved by Sachs and Wolfe (§ sec. II), confirmed independently by Lukash |^], 
Chibisov and Mukhanov on the basis of Hamilton formalism in Field-Shepley variables Q]. The result holds 
in the gauge-invariant descriptions as well as in the original Lifshitz formalism 

The situation changes significantly in the universe with the negative space curvature {K = — 1), where the 
acoustic wave propagates as a scalar field with the mass m = — itT.Q The dispersion relation is nonlinear and, in 
consequence, some critical frequency Wcr appear. Below ojc the wave propagation is forbidden. The dispersion 
of the acoustic wave on the curvature affects the microwave background fluctuations and in principle can be 
used to measure the cosmological density parameter 12 . 

In this letter we show that also the cosmological constant H is a source of dispersive phenomena. The 
interest in H-cosmology grows, owing to speculations that cosmological constant may take substantial values 
in the present epoch |^. We investigate the flat [K = 0) universe with radiation {p = e/3) and positive 
cosmological constant (T > 0). We And exact solutions for gauge-invariant perturbation equations |^--[|^ and 
discuss their properties to show, that in the presence of positive A the dispersion of the acoustic field has an 
anomalous character. 


2 Homogeneous background 


In the flat radiation-dominated universe with a positive cosmological constant, Friedman equations are satisfied 
by the scale factor a{t): 


a{t) = {Mj \J sinh(2\/yl/3t) = {Mj Vsinhr (I) 

where M. is the constant of motion AA = e a^, e means the energy density and the dimensionless time parameter r 
is related to the metric time as t = 2^ AjZt. To simplify notation in the forthcoming formulae we also introduce 
the constant /3 defined as /3 = •\/3/[16A^yl]^/^. The conformal time 77 


r] = F 


^arccos 


I — sinh T 

1 ^ 

I -|- sinh T 



( 2 ) 


is finite ry G [ 0 , 2iF(2“^/^)/3] = [ 0 , i 7 r“^/^[r(i)]^/ 3 ], where F{p, tu) stands for the elliptic integral of the first 
kind and K{m) is the complete elliptic integral of the first kind ||^. The scale factor 0(77) expresses by the 
Jacobi elliptic functions sn{u,m), cn{u,m) of the ratio rj/p. 


sn(:2 

air,) = 2Py^ (3) 

The scale factor 0 ( 77 ) grows monotonically from 0 to infinity. In the neighborhood of the initial singularity the 
radiation is dynamically dominant and the effect of the cosmological constant can be omitted. The opposite 
boundary of the conformal time interval (77 = i 7 r“^/^[F(;|)]^/ 3 ) correspond to the de Sitter epoch where A 
dominates and the radiation plays a marginal role. 

^Gauge-invariant perturbation equation takes the form of eq. (5.4) of with the minimal coupling (g = 0). The same is also 
the form of eq. (21) of and eq. (4.6) of 
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3 The perturbation equation 

By variating Raychaudhuri and the continuity equations one obtains the propagation equation for the density 
contrast 5e/e. While 5e/e is measured on the flow-orthogonal hypersurfaces |Q the equation reads 

— i coth^(T)X(r, ^ coth(T)X'(r, x) -|- X" {t) = v‘^0^ csch(r)V^X(T, x). (4) 

denotes the Laplace operator in the conformal space, consequently its eigenvalue —k'^ is related to the time 
invariant “commoving” wave vector k with k = |k|. The same equation can be obtained in other gauge-invariant 
theories ^-[|ll|, under necessary redefinition of the perturbation variable^ (see |^). This is convenient 

to express the equation in terms of the conformal time 

7^ -1 ^ ) X{r],x) - v^\/^X{7],x) + X''{r],x) =0 (5) 

and search for its solutions in the form of the Fourier integral 


X{r],x) 


j ylk Uk(? 7 , x) -I-c.c. 


( 6 ) 


Prime stands for the conformal time derivative, Fourier coefficients Ak are arbitrary functions of the wave 
number k. Modes Uk(r 7 ,x) satisfy both, the Helmholtz equation 


V^Uk(??,x) = -fc^Uk(ry,x) 


(7) 


and the time equation (i , and divide into two classes distinguished by the critical value k = = (4Af 

Below, we employ the dimensionless wave vector |k| = k = \/2v(3k = In this notation both classes of 

solutions are defined by inequalities k > 1 and k < 1. They will be discussed separately. 


4 Subcritical perturbations: solutions for «: > 1 

The short-scale solutions (k > 1) propagate as waves of variable amplitude and variable frequency. Elementary 
solutions Uti{r],x) take the form 

sn([2v^/3]-i[r(i)]2, -^)J2 + {k^ - l)sn2(^,^) 

y - -exp[i 6>(77 ,x)] (8) 

2 + _ 1) sn2([20F/3]-i[T(i)]2, 

where the phase 0 (? 7 , x) is given by 



0{r],x) 


K • X 

y/2f3v 


Ky/lA 

V2{k^ 


1 

~) 


V 






(9) 


Function 7j(v5, m) is the elliptic integral of the third kind and am(u, m) is the amplitude elliptic function [T^ . 
The solution is singular at both ends of the conformal time interval. In the vicinity of the initial singularity 

^Although we restrict ourselves to v = I/'s/s we explicitly keep v in perturbation equations to make clear their structure and 
correspondence to other formalisms. 
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{t] = 0) the radiation plays the dominant role and the solutions Uk(? 7 ,x) approach appropriate modes for the 
universe with vanishing A In the late time regime t] ~ universe realizes the exponential 

growth under the dominant effect of A. Singularities in solutions Uk(? 7 , x) in the r/ —> limit is known 

as the instability of de Sitter universe Q . 

The local frequency^ evaluated as the time derivative of the phase 0{r], x) 


d0{T],K) KVK‘‘-lsn 2 (|,^) 

dv " V2/3(2 + (K 2 _l)sn 2 (^,^))- 


( 10 ) 


is a time-dependent quantity. The interval rj G [? 7 i, 772 ] contained between two extrema of the second derivative 

d‘^uj{r},K.) 

dr]"^ 


7]i = F ^arcsin 



V2 = \7^ - F (3 


form the oscillatory stage, where the perturbation propagates as the acoustic wave. w(? 7 , k) maintains an 
approximately constant value there. The quantity 


CJb(K) 


kV — 1 
V2/3(At2 + l)' 


( 11 ) 


may be understood as the basic frequency, while its the modulation oJmodiv^ ^) is given by 


^^mod(?7>K) 


V^kVk^ - 1(-1 + sn 2 (^, ^)) 

/3(l + «2)(2 + (^2_l)sn2(|^^))- 


( 12 ) 


Beyond this interval [r]i, 772 ] the division into basic frequency and the modulation loses its physical meaning. 
While 77 —!■ 0 or ?7 ^ 11^® frequency formally defined in ( 0 ) tends to zero and the amplitude grows 

indefinitely. In fact, the concept of oscillation, breaks down in these limits. 

Formula (^^ plays the role of time-dependent dispersion relation, and ( |Tl| ) is its analog for the basic 
frequency. Both relations are nonlinear. The number k = 1 {k = ^273/3 = forms the critical value 

of the wave number, and naturally defines the critical wave-length Acr = Below this scale the wave 

propagation is forbidden. While k the frequency tends to zero, therefore, perturbations longer than the 
critical size do not form travelling waves. In contrast to dispersion on the space curvature, the dispersion on yl 
manifests its anomalous character: critical behaviour concerns the scale of the perturbation (the wave number 
k), not the frequencyQw. The anomaly is also clear from the phase and the group velocities behaviour 


;(77,K,/3) = 


Vi{v,K,P) = 
907(77, k ) 


}{r],K) 




d k 


k V3[2+(«2-I)sn2(f,^)]’ 

[6k 4 _ 2 -h (k 2 _ 1)(«;4 _ 2^2 _ ■^)] sn2(|, 

V3V)^^[2 + (7«2-l)sn2(|,^)]2 ■ 


(13) 


(14) 


The phase velocity 77 f ( 77 , k, (3) decreases with increasing wave-length and, in consequence, the group velocity 
Vg(r],K,f3) is always greater than the phase velocity V{{r], k, (3). In the k 1 limit the phase velocity tends to 
zero, while the group velocity formally expressed as ([l^ ) grows indefinitely. In this regime the concept of the 
wave packet loses its sense.| 

^see ligj 

^The space curvature defines the minimal frequency like the plasma frequency. 

®Taylor aeries for ui cannot be cut at linear terms. Higher order terms are responsible for the diffusion of the wave packets 0- 
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5 Supercritical perturbations: solutions for k, <1 


The large-scale solutions {n < 1) do not propagate as travelling waves. The space of solutions consists of the 
combinations of the two modes exp(z«; • x) and exp(i«; • x), where both and 




sn([2V^/3]-i[r(i)]2, ^2-(l-«2)sn2(|,^) 

^2-(l-«2)sn2([20F/3]-i[r(i)P,^) sn(f, 


exp 


K^/l — 

V2(l - k2) 








(15) 


u 


( 2 ) 


iv) 



1 


dri'. 


(16) 


are real functions. In the vicinity of the initial singularity behave as the decaying mode, while plays 
the role of the growing one. In this particular sense the behaviour of the large-scale perturbations mimics that 
of modes in the pressureless matter. Later on the solutions evolve differently. Both modes keep approximately 
constant amplitude in the epoch when the k > 1 modes oscillate. Finally, in the de Sitter regime both modes 
and blow up exponentially. 

The two limit cases k = 0 and k = I express in simple formulae 




sn([20F/3]-i[T(i)]2,^) 


2 — sn" 


(d’ 72 ) 


2-sn2([20F/3]-i[T(i)]2,^) sn 
sn([20F/3]-i[T(i)]2,^) 
^2-sn2([20F/3]-i[T(i)]2,^) 


(d’ V 2 ) 


2 — sn" 


(^’ V 2 ) 


v^) 


V 1 


2^^ am 


1 


P’V2j’V2 


V 


P 


Tj 1 

-"'(is'Ti 


for K = 0 and 


uW(? 7 ) = sn 
^k\v) = sn 


(^[2V^/3]-i[r(i)]2, 


V 2 J sn(|,^)’ 


1 


V2J sn(f,^) 


V 



V J_ 

P’ V2 


1 

’72;j 


(17) 


(18) 


(19) 

( 20 ) 


for K = 1. 

The dynamics of large-scale perturbations does not confirm their exceptional role in the structure formation 
process: their amplitudes are of the same range as those of oscillating modes. However, since these perturbations 
do not propagate, the density and expansion excesses are correlated to each other throughout the space for 
decaying modes, and anti- correlated for the growing ones. Present observations of the microwave background, 
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which are limited to the single physical quantity (temperature), are not complete enough to distinguish between 
traveling waves and the perturbations of other types. Yet it is possible that we will get more complete data from 
the last scattering surface in future. Therefore, some estimations of linear diameters and appropriate angular 
scales on the sky are worth mentioning. 

We assume that the universe after decoupling evolves as filled with presureless matter. Then the lapse of 
conformal time between the emission and the observation moments is given by 


Ari{z) = \/i/A 


fo 

LmJ 


nl/3 


(1 + Z)2i^l 


- -• -■ -(1 + z)'^ — 
3’ 2’ 3’ ^ ’a 


— 2F1 


1 1.4 eo 
3’ 2’ 3’ A 


( 21 ) 


where z means the redshift of the last scattering surface, m = e is the constant of motion for the late universe, 
subscript 0 refers to the present time, and 2 ^"i[a, b; c; x] stands for the hypergeometric function |2|]. Sky equator 
covers n regions of the critical length diameter where 




27rAry 


= Ve 


eo 




1/4 


(1 + z) 2 F 1 


i -■ -■ _L ,'|3£o 

3’2’3’ ^ A 


— 2 F 1 


1 14^ 
3’ 2’ 3’ 


£0 

'A 


( 22 ) 


For cq!A = 0.2 and z = 100 the number n is about n = 0.8, therefore, supercritical perturbations may contribute 
to the dipole anisotropy of the MBR. If the radiation plays its dynamical role after decoupling the number n 
can be greater, but always less than n = 2. 


6 Conclusions 

Like the space curvature, the cosmological constant causes dispersion of the acoustic field in the early universe. 
In both cases the dispersive phenomena are induced by the space-time geometry and depend neither on the 
initial conditions for acoustic field nor on the character of cosmological epochs prior to the radiational era. 
Dispersion on the cosmological constant is a measurable effect in the sense that appropriate critical length scale 
is comparable with the region of the last scattering surface we observe today. Nevertheless, to detect dispersion 
one needs more complete data from the last scattering epoch, in particular, the second dynamically independent 
quantity should be known (for instance the velocity field). In contrast to dispersion on the space curvature, 
the dispersion on A is anomalous. The two phenomena differ enough to uniquely distinguish their provenience. 
Both might be employed to measure values of A and 12 respectively. 


Acknowledgements 

This work was partially supported by State Committee for Scientific Research, project No 2 P03D 014 17. 

References 

[ 1 ] R.K. Sachs and A.M. Wolfe, Astrophys. J. 147, 73 (1967). 

[2] V.N. Lukash, Sov. Phys. JETP 79, 1601 (1980). 

[3] G.V. Chibisov and V.F. Mukhanov, Mon. Not. R. Astron. Soc. 200, 535 (1982). 

[4] G.B. Field and L.C. Shepley, Astroph. Space Sci. 1, 309 (1968). 

[5] Z.A. Golda and A. Woszczyna, Class. Quantum Grav. 18, 543 (2001). 


6 
















[6] Z.A. Golda and A. Woszczyna, J. Math. Phys. 42, 456 (2001). 

[7] N.D. Birrell and P.C.W. Davies, Quantum Fields in Curved Space (Cambridge University Press, Cambridge, 
1982). 

[8] N.A. Bahcall, J.P. Ostriker, S. Perlmutter and P.J. Steinhardt, Science 284, 1481 (1999). 

[9] J.M Bardeen, Phys. Rev. D 22, 1882 (1980). 

[10] R. Brandenberger, R. Kahn and W.H. Press, Phys. Rev. D 28, 1809 (1983). 

[11] D.W. Olson, Phys. Rev. D 14, 327 (1976). 

[12] D.H. Lyth and M. Mukherjee, Phys. Rev. D 38, 485 (1988). 

[13] A. Woszczyna and A. Kulak, Class. Quantum Crav. 6, 1665 (1989). 

[14] C.F.R. Ellis and M. Bruni, Phys. Rev. D 40, 1804 (1989). 

[15] C.F.R. Ellis, M. Bruni and J. Hwang, Phys. Rev. D 42, 1035 (1990). 

[16] J.C. Jackson, Mon. Not. R. Astron. Soc. 264, 729 (1993). 

[17] P.F. Byrd and M.D. Friedman, Handbook of Elliptic Integrals for Engineers and Physicists (Springer-Verlag, 
Berlin, 1954). 

[18] W. Zimdahl, Mon. Not. R. Astron. Soc. 257, 445 (1992). 

[19] C.B. Whitham, Linear and Nonlinear Waves (John Wiley & Sons, New York, 1974). 

[20] Sz. Szczeniowski, Experimental Physics (Part I Acoustics and Machanics), (in Polish) (PWN, Warszawa, 
1972). 

[21] W.L. Cinzburg, Rasprostranenie ehlektromagnitnykh voln v plazme, (in Russian) (CIF-ML, Moscow, 1960). 

[22] I.S. Cradshteyn, I.M. Ryzhik, Table of Integrals, Series, and Products 5th edn (Academic, New York, 
1994). 


7 



